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Interacting electrons with a square Fermi surface is investigated from a bosonic point of view 
taking into account electron scattering between all faces of the square. Fermion operators are 
classified according to their dimensions and the stability of the boson fixed-point is investigated. 
In particular we find, in contrast to previous studies, that the square Fermi surface is unstable to 
doping in the case of no spin gap and microscopic Hubbard interactions. 



PACS numbers: 71.10 Fd 



I. INTRODUCTION 



Much effort has been invested in solving the two di- 
mensional one-band Hubbard model which describes re- 
pulsive electrons hopping on a square lattice. The poten- 
tial payoff is a possible explanation of high-Tc supercon- 
ductivity. Yet no definite solution exists and it might be 
more valuable to shift focus from the microscopic Hub- 
bard model to more general classes of effective models. 
Such a strategy is well formulated using the renormal- 
ization group approach. The essence of this approach 
is that low-energy properties of any microscopic Hamil- 
tonian become identical, at low enough energies, to the 
same properties of an effective Hamiltonian describing 
a fixed point of the renormalization group flow. Insta- 
bilities and phase transitions can be understood by the 
presence of relevant operators that cause flows away from 
this fixed point. 

The best studied fixed point for 2D electrons is that 
of free electrons with a Fermi surface. There are only a 
few instabilities associated with this fixed point*, among 
them the BCS superconducting instability. Another fixed 
point for 2D electrons can be arrived at by considering 
electrons on a square Fermi surface as a collection of 
electrons on coupled chains 2 . This fixed-point Hamil- 
tonian consists of the chain electron density operators as 
well as forward-scattering terms and is purely bosonic. 
In the context of 2D electrons this boson Hamiltonian 
has not been extensively studied except for the works 2 ** 
where a simplifying assumption was employed: No inter- 
actions between density operators on adjacent faces of 
the square. 

In this article we extend the analysis of the boson 
Hamiltonian to also include adjacent face couplings in 
such a way that the most general interactions respecting 
the square symmetry are taken into account. We diag- 
onalize explicitly the boson Hamiltonian, and classify a 
family of fermion operators with respect to this fixed- 
point. As will be shown, a particularly important conse- 
quence of adjacent face interactions is that the stability 
of the square Fermi surface is reduced. 

It should be mentioned that the same boson Hamilto- 
nian as considered here, with other values for the param- 
eters, appears in the studies of the so called sliding Lut- 



tinger liquids*, both in the context of modeling striped 
phases of cuprates* and the quantum Hall effect. 



II. THE SQUARE FERMI SURFACE MAPPING 

The construction used to map fermions on a square 
Fermi surface to chains was explained in details in Ref^. 
In the following paragraphs we will outline the essentials 
of this mapping. The construction is similar in spirit 
to what is done in one dimension. There the real-space 
fermion operator is written as a sum of two terms, each 
dominated by momentum components around one of the 
Fermi points. On the square Fermi surface, the fermion 
operator is decomposed into a sum over the four faces of 
the square 
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where ip p has momentum components coming from face 
p. This face field can be written in terms of momen- 
tum components, which have a natural decomposition 
into momentum components parallel and perpendicular 
to the face normal. Identifying the momentum compo- 
nents having identical perpendicular momentum with the 
momentum components of a one dimensional chain close 
to one of its Fermi points, one can write the face field 
as a sum over real-space chains oriented along the face 
normal 

M*>v) = \hfY,9{i-i')i>i'Ax) (2) 
V / i> 

where kf — 7r/(v2a), and 

N/2 

2N 



n y = -N/2 

is the real space structure function coming from restrict- 
ing the perpendicular momentum coordinate which la- 
bels each chain, to lie within [—kf, kf]. We have rotated 
the coordinate system such that the integer label I is re- 
lated to the real space coordinate y parallel to face 1, 
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FIG. 1: The square Fermi surface 



y/2y = al. where a is the lattice spacing. Note that the 
square Fermi surface is used here as a bookkeeping tool 
only in order to keep track of relevant states. Thus the 
real Fermi surface of electrons does not necessarily need 
to be square. 

The Hamiltonian is obtained by first considering the 
kinetic energy term of free fermions with a square Fermi 
surface. Linearizing the energy spectrum and using the 
mapping Eq. @ one obtains a bilinear term of fermion 
operators on in general different chains, but belonging 
to the same face. Only those terms where the fermion 
operators belong to the same chain is retained in the 
fixed-point Hamiltonian, thus the fixed-point Hamilto- 
nian describes effectively a collection of ID chains. This 
is bosonized in the standard way giving a Hamiltonian 
quadratic in density operators for each chain. The effec- 
tive velocity vq associated is an average over the fermi 
velocity over the face. The rest of the terms involving 
different chain operators are treated as perturbations to 
the boson fixed-point Hamiltonian. There are also other 
terms corresponding to fermion-fcrmion interactions that 
can be written as a product of two chain density opera- 
tors. Using the mapping Eq. @ for four- fermion interac- 
tions, pairwise chain contributions corresponding to pairs 
of density operators are also included in the fixed-point 
Hamiltonian, giving the fixed-point Hamiltonian proper- 
ties of the interacting fermion system. 

The most general such Hamiltonian invariant under 
rotations quadratic in density operators is 



H 
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Ppk± (^p )A pp i (k p , k p i )p p _ k ± ( fcjj ) (4) 



a—s,c 



where the labels p and p' denotes the faces of the square 
and takes values from 1 to 4. The momentum k p is the 
momentum parallel to the face normal, that is along the 
chain and k p is the momentum along the Fermi surface, 
perpendicular to the chains. We have divided the vari- 
ables into charge and spin sectors: v / 2p c = P] + p[ and 
= PT - Pi ■ Tne matrix entries of A are restricted 
by the rotational invariance. We have An = A + + A~ , 
Aij = A + — A~ for i and j labeling opposite phases, and 
Aij = V for i and j labeling adjacent faces. 

The contributions from the kinetic energy to these cou- 



pling constants are 
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Mdn = 

being the same in both the spin and the charge sector. 
We will treat the Hamiltonian for generic couplings, but 
it is illustrative to write down the explicit expressions 
in the case when the interaction term is the Hubbard U 
term 



U 



Then the charge sector parameters take the form 
Al = ^ + ^f(k-) 



(5) 



At = 



TTVQ 



(6) 



V c = ^/(fc 11 )/^) 



where f(k) is the Fourier transform of g 2 given by 

/(*) = xyco^ = 5 (i - a) 



The spin sector parameters are obtained by letting U 
-U. 



III. DIAGONALIZATION 

In order to diagonalize the Hamiltonian Eq. Q we 
write it down explicitly. It has the same form in the two 
(a = c, s) sectors 



2NL 



A 



'hi 



{pig v {q x ) + vp 2qy (q x j) 



(8) 



+ B t : I (Pzq* (Qv) + v Pii* (Qv)) I 2 : 

+ V q,q y ■ [(Plqy feO + P2q y (&)) 

x (P3q* (Qy) + P4q* (?y)) T + h.c. :| 

where we have used B to denote the function A when 
its argument is q x , and V qxq = V r gy9x . For notational 
ease we have omitted the superscript a on the coupling 
constants. This is the most general Hamiltonian bilinear 
in density operators consistent with the symmetries of 
the square (A^ can also depend on q x generally, but this 
is omitted here for simplicity). 

It is convenient to introduce the combinations 



Pi ± P2 



P3 ± Pi 



J^(e(g x )±e(- qx ))(all«±a-_%) 



L \%\ 
2tt 



m v )±9(-q v ))[b% r ±bZt 



(9) 
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where 9 is the step-function and a and b are commuting 
boson annihilation operators. In order to diagonalize the 
Hamiltonian we express the a and b operators in terms 
of new operators a and /3 as follows 

a%>±aZ% = e^cos 7 (atf ±oq:) 

+e±* sin 7 (/^± (10) 

blT±b-j q ; = e ^'cos 7 (e±/q;) 

-e**' 'sin 7 (a£"±a:£) (11) 

where j,0,0', (j),(j)' are real numbers. It is con- 
venient to shorten the notation somewhat writing 
A ± = Af v L\q x \ /(2n), B ± = B+ L\q y \/(2w) and V = 

Vq,,q v L V\Q^\\QV\/( 2 ^)- 

The a's and /3's commute with each other. The preser- 
vation of the commutation relations of the a's and b's 
imply the condition 







(12) 



By inserting the above expression into the Hamiltonian 
one finds a term proportional to (o^ — a)(ft — [3). The 
coefficient of this term can be set equal to zero by setting 



o 2(0+0') _ 



B- 



(13) 



Similarly there is a (cJ +a)(ft +(3) term which coefficient 
can be set to zero by choosing the angle 7 to fulfill 



tan(2 7 ) = 



2|vivnr 

B+\B~\-A+\A- 



(14) 



There are also a) a) terms (similarly for /3), getting rid 
of these fixes the angles 



o 20 



26' 



e 20' = 



/ (^-) 2 (l-tan 2 7 ) 
B+B- -A+A- tan 2 7 ' 



-20 



/ (g-) 2 (l-tan 2 7 ) 
A+A- -B+B- tan 2 7 ' 



26 _ B -2<t>' 

A- ■ 

Having fixed the angles the Hamiltonian is diagonalized 
and reads 



(15) 



where 



K = g a + (A+A- cos 2 7 + B+B~ sin 2 7 ) /g a , (16) 



and the parameter g a = <J A+A i_taJ 7 ta " 27 - h P is g ot " 
ten from the expression for h a by interchanging *4s and 
Bs. 

The asymptotic behavior of correlation functions are 
governed by small momentum values. For \q x \ -C \q y \ we 
find 
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h fi = 




B o B o ■ 



The behavior for \q y \ <C \q x \ can be found by interchang- 
ing A <-» f?,^ <-» g y and primed <-» unprimed quantities. 



IV. OPERATOR DIMENSIONS 

We will now consider the stability of the square Fermi 
surface fixed point to products of fermion chain operators 
on equal and opposite faces: 

° = ^Vinh^P^-ih ■ ■ ■ ^PmtmIm ( 17 ) 

all acting at the same x and r, and where D = +1 
means a creation operator and — means an annihila- 
tion operator, p = +1,-1 corresponds to face 1 and 
2 respectively and a = ±1 denotes the spin direc- 
tion. To write the expression for the operator dimen- 
sion of the above operator in a compact form, we in- 
troduce four integer parameters, q$j = ^2m—iSi,i m D m , 
q<t> c l = J2m=ih,i m PmD m , qej = J2m=ihjm a mD m , and 
*M = Em=i 5i.i m p m a m D m . The dimension d defined by 
(O(x)O(0)) = l/x 2d is then 



d = 



1 1 

427 



J-7T Z „, 



qejqej' (18) 



+K ck ±q l pjq< l , c i> + K^qejqgj, + K sk j_q^jq^ji) 

where I, V are integers labeling the chains such that neigh- 
boring chains differ by unity. The sum runs over all the 
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TABLE I: Table of operator dimensions. The dimension is 
d = J ^d(k) where d(k) is gotten from the table. The 
abbreviations (+) = 1 + cos(fcZ) and (— ) = 1 — cos(fcZ). An 
asterisk * indicates that the operator does not correspond to 
a term in the microscopic Hubbard model. 



# Operator 


d(fc) 




i ; t / 

1 


i^ + a-^r + ka. + a^X-) 




(Ac + A- 1 )^) 




(A s + 






A- x ( 


-) + AT 1 (-) 


r /,t / T „ / / 


A- x ( 


-) + K s {+) 


6 ?i'J-»i/'l- i thut ibin i 


A" x ( 


-) + K s (~) 




K-'i 


-) + K s (+) 




K c (- 


) + A s (-) 




K c (- 


) + K s (+) 


10 ^I^^V'ijT'^a 


K c {- 


) + K- 1 (-) 




K~\ 


-) + K c (+) 


12* VLt^Iji^t^-t 


A c (+ 


) + K s (+) 




A c (+ 


) + K- 1 (-) 


14 VutWil^aiT^yi 


A c (+ 


) + *.(-) 



TABLE II: Dimensions for operators occurring in various 
response functions. The dimension is d — ^d(k) where 
d(k) is gotten from the table. The abbreviations are as in 
Table □ 



type 


Operator 


d(fc) 












U K cH+) 


f A c (-) 




+ A S (+ 










+■ K-\+) 


+ A S (- 


CDW 






+■ A c (+) 




+ A S (+ 


SDW 




\{K-\-) 


+■ A c (+) 


f A7 1 (+) 


+ A S (- 






\{K-\+) 


f A c (+) 


f A s - 1 (+) 


+ A S (+ 








f A c (+) 




+ K s (- 




V'LtV'ijI. 




fATc(-) 


f A s - 1 (+) 


+ A S (+ 



chains present in the operator. The Luttinger parameter 
which depends on the perpendicular momentum is given 
by 



KuA 



[cos 2 ( 7 )e 2e +sin 2 ( 7 )e- 2 ^ 
1 




fel!->0 



(19) 



This expression for the Luttinger parameter was also 
found in the context of crossed sliding Luttinger liquids 7 . 

Specializing to the Hubbard Hamiltonian Eqs. 10, the 
Luttinger parameter in the charge sector K c can be writ- 
ten 



K c (k x ) 



/ 2 (fc J 



(20) 



where u = and / 2 (0) = 1/4. the similar expression 
for K s differs by having u — > — u. The change introduced 
by the adjacent face couplings is thus to renormalize the 
dimensionless coupling u — > u/ (1 + uf 2 (Q)). While this is 
a small effect for small u, it is very important for large u 
in that it isn't possible to lower the value of K c arbitrarily 
much by increasing u. In fact the lowest possible value 
of K c is l/\/2 which is reached for k 1 - — at u — > oo. 
For k 1 - — > 7T the value of K c approaches 1 for all values 
of u. 

Table [I] shows the operator dimensions of all spin- 
conserving two and four-fermion interactions on same 
and opposite faces. Operators that reduce to pure den- 
sity operators are not written down. A dimension smaller 
than 2 implies that the operator is relevant. Similarly Ta- 
ble^ shows the operator dimensions of operators occur- 
ring in various response functions including singlet (0 SS ) 
and triplet ( Ots) superconductivity as well as charge and 
spin density waves (CDW and SDW respectively). 



V. RELEVANT OPERATORS 

While K s is different from unity for the Hubbard in- 
teraction, the SU(2) symmetry in spin space requires K s 
to flow to 1 as long as the spin-correlators decrease al- 
gebraically. However if the system develops a spin gap 
the correlators decay no longer algebraically and the spin 
sector is not part of the low energy physics. Thus we con- 
sider these two cases separately. 



a) No spin gap 

For the square Fermi surface to be stable upon doping, 
the operator V'jijV'jiT which comes from the chemical po- 
tential should in particular be irrelevant, that is its di- 
mension should be greater than 2. Using the expression 
Eq. (|18fl inserting K s = 1 we find that the condition for 
irrelevancy is 

1 f V rlh 

4 / ^(^( fc )+^ 1 (fc) + 2)(l-cos(^-j)))>2. 

(21) 

For momentum independent K c it is easy to see that 
the chemical potential term is relevant for 3 — 2^/2 < 
K c < 3 + 2y/2. Taking into account the weak momentum 
dependency given by f 2 (k) and focusing on the near- 
est neighbor chains such that i — j = 1 we find numer- 
ically that the chemical potential term is irrelevant for 
u/(l + u/4) > 1378.76. This inequality is clearly not 
fulfilled for any positive values of U, thus we conclude 
that the square Fermi surface is unstable to Hubbard- 
like interactions upon doping. This is at odds with the 
conclusions arrived at in Refsp2*2> where it was found that 
a sufficiently large value of U would make the chemical 
potential term irrelevant. However there adjacent face 
couplings were not taken into account. 
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FIG. function of the anisotropy parameter ^y a for 

two values of 70. For u > u c , the chemical potential term is 
irrelevant and the square fermi surface is stable upon doping. 



It is interesting to investigate what the anisotropics be- 
tween the different face scattering rates must be in order 
for the square to remain stable. Introducing anisotropy 
parameters describing the ratio of opposite to same face 
interactions 7 = t/12/^11, and adjacent to same face 
interactions 7 a = U13/U11 we can write 



K. 



1 



ck 1 



P{k^)} [1 + 



,1+7° 



[1 



,1+7, 



-/ 2 (fc^)] [1 + 



,l±7o- 



f7^/ 2 (fc ± ) 

where we have used / 2 (0) = 1/4. For small u, K ck ± is 
close to unity and the inequality (|21|l is violated. Thus 
this stability criterion sets a lower limit u c on u. For 
big it, K c (x \/l — 7 , thus 70 should be close to 1 in 
order to get a small K c which can saturate the inequality. 
How close can be found by setting u — > 00. It follows 
that 7„ > 2V2/3 m 0.9428. We have used the fact that 
7 Q generally increases K c , thus the lowest U c occurs for 
7 a = 0. In Fig[21wc is plotted as a function of 7 a for two 
values of 70. Note the extreme values of u needed for 
the square to remain stable regardless of the anisotropy. 
It is unclear how such extreme values can emerge in real 
systems of strongly correlated electrons. 

The chemical potential term is not the only rele- 
vant operator when K c behaves according to Eq. H2()(l 
taken from the microscopic Hubbard model. As seen 
from Table [I] other relevant operators are the back- 
scattering interactions (8)-(10) and umklapp processes 
(12)-(14). There are also interactions that affect the 
spin sector only: Operators (3) and (7) can be inter- 
preted to arise from interactions between spin variables 
on same/opposite faces and different chains. In the 
isotropic case K s = 1 they are both marginal operators. 



b) Spin gap 

When the system develops a spin gap, spin correlators 
do not decay as power-laws and the value of K s is not 



fixed to unity anymore. Formally this situation can be 
covered by setting K s — > or 00 dependent on which 
field is pinned down. In either of these cases, because 
the chemical potential term involves both K s and K~ x 
it will automatically be irrelevant^. The physical reason 
for this is that removing a particle from a chain involves 
breaking the spin gap, and is therefore not part of the 
low-energy physics. In the case where K s — > the opera- 
tors (5), (6) corresponding to singlet pair- hopping, (8), (9) 
to different chain back-scattering, and (12), (14) umklapp 
interactions are all potentially relevant dependent on the 
expression for K c . On the other hand when K s — > 00 the 
operators (4) corresponding to triplet pair-hopping, (10) 
to different chain back-scattering and (13) to umklapp 
scattering are potentially relevant. The value of K s got- 
ten from the microscopic Hubbard model indicates that if 
a spin gap forms it should be of the form K s — > 00. Using 
this together with Eq. (J20J operators (4), (10) and (13) 
are all found to be relevant when considering nearest- 
neighbor chains. 

The operators (2) corresponding to same-direction pair 
hopping and (11) umklapp scattering do not depend on 
K s . With K c following Eq. l|2Tj|) operator (2) is found to 
be irrelevant while (11) is relevant for nearest-neighbor 
chains. When taking also inot account operators that de- 
scribe interactions between chains at bigger separations 
(|i — j\ > 1) it becomes clear that there are a number 
of potential instabilities and a careful renormalization 
group treatment which depends on the details of the mi- 
croscopic Hamiltonian is needed to determine the domi- 
nant instability. 



VI. CONCLUSION 



The basic insight gotten from taking adjacent face in- 
teractions into account in the boson Hamiltonian is that 
it is not possible to make K c deviate significantly from 
unity even for large values of the Hubbard-U interaction. 
Thus the ultimate fate of the system depends to a large 
extent on what happens in the spin sector. If no spin gap 
develops, the square Fermi surface is unstable upon dop- 
ing unless there is a significant anisotropy between adja- 
cent and same face scattering processes and the value of 
U is very large. On the other hand if a spin gap develops 
the chemical potential term is irrelevant. However in that 
case there are also a number of other relevant operators 
which makes the dominant instability dependent on the 
details of the microscopic Hamiltonian. 
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